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Abstract — Phase noise and frequency offsets are due to their 
time-variant behavior one of the most limiting disturbances in 
practical OFDM designs and therefore intensively studied by 
many authors. In this paper we present a generalized frame- 
work for the prediction of uncoded system performance in the 
presence of time- variant distortions including the transmitter and 
receiver pulse shapes as well as the channel. Therefore, unlike 
existing studies, our approach can be employed for more general 
multicarrier schemes. To show the usefulness of our approach, we 
apply the results to OFDM in the context of frequency offset and 
Wiener phase noise, yielding improved bounds on the uncoded 
performance. In particular, we obtain exact formulas for the 
averaged performance in AWGN and time-invariant multipath 
channels. 

Index Terms — Multicarrier transmission, OFDM, Gabor the- 
ory, frequency offset, phase noise 



I. Introduction 

Multicarrier (MC) transmission is a promising concept 
for future mobile communications. In particular, the popular 
OFDM scheme has already been implemented in several stan- 
dards such as digital video and audio broadcasting (DAB and 
DVB-T), WLAN (IEEE 802. 11a) and is proposed for next gen- 
eration mobile communication networks. The main advantage 
of OFDM is its easy implementation and simple equalization 
which is based on the orthogonality of the subcarriers in time- 
invariant channels. However, time variances due to Doppler 
shifts, carrier frequency mismatching or phase noise destroy 
this orthogonality and causes intercarrier interference leading 
to serious performance degradation in OFDM transmission [1]. 

Most authors [2], [3], [4], [5], [6] study offsets and phase 
noise in the context of OFDM. Unfortunately their approaches 
are difficult to apply to many recently proposed MC schemes, 
which are optimized in terms of pulse adaption and bandwidth 
efficiency [7], [8], [9], [10]. Even in the OFDM case exact 
results are not known. However, exact results are crucial for 
next generation systems, that are expected to operate at high 
mobility as well as at high data rates. Moreover, to overcome 
the disadvantages of low-cost hardware at the mobile side 
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normally some kind of tracking is used to provide a minimum 
signal quality. On the other hand tracking increases receiver 
complexity. Hence in performance evaluation this has to be 
taken into account which is provided by this contribution. 
Another goal of the paper is to establish a new approach 
for general MC schemes without the restriction to "textbook" 
OFDM. The main idea is, that for many practical schemes 
the average total power received on a single subcarrier is 
preserved, giving bounds on the expected interference power. 
Hence, interference analysis reduces to the investigation of 
separated subcarriers. 

The paper is organized as follows. First the cyclic prefix 
based OFDM transmission model which is relevant for most 
applications is introduced. Next the transition to generalized 
MC transmission is performed together with some introductory 
notes on the underlying Gabor theory. For this system model 
we study in Sec. Ill the effective mapping that results from 
time-invariant channel and a linear distortion. The evaluations 
lead to a theorem on the resulting interference. In Sec. IV and 
Sec.V the results are applied to time-frequency offsets and 
phase noise. 

II. System Description 

A. OFDM signaling 

The cyclic prefix based OFDM (cp-OFDM) baseband trans- 
mit signal is 

T 

s(t)= Yl x^e^/^jit-nfrT+T^)) 

(mn)£l 

where i is the imaginary unit and 7(-) is the rectangular pulse: 

1 {t) = -===X[-T cp ,T u ]{t + t ). (1) 
y/ -*u < -L cp 

The function X[-t cp ,t u ] is the characteristic function of the 
interval [— T cp , T u ], where T u denotes the length of the useful 
part of the signal and T cp the length of the cyclic prefix. 
Without loss of generality we set the time origin to to = 0, but 
note that this has influence on several common phase errors 
given in this paper. The subcarrier spacing is F = 1/T U 
and x mn are the complex data symbols at time instant n 
and subcarrier index m. The indices (mn) range over the 
doubly-countable index set I, referring to the data burst to 



be transmitted. Note that in practice only a finite number of 
subcarriers is considered. However for theoretical reasons it 
is beneficial to consider an infinite set. We will denote the 
synchronization defects by the random linear operator S, the 
linear time-invariant channel by Ti. and the additive white 
Gaussian noise process (AWGN) by n(t). The received signal 
is then 

r(t) = {SUs)(t) + (Sn)(t) 
= (S(s * h))(t) + (Sn)(t) 

with h being a realization of the (causal) channel impulse 
response of finite maximum delay spread r^. The standard 
OFDM receiver estimates the complex symbol as 

x k i = J e-^ kt ^g(t-l(T u + T cp ))r(t)dt 

using the rectangular pulse g{t) = -±=X[o,T u ]{t + *o) which 
removes the cyclic prefix. If we assume that the receiver 
has perfect channel knowledge (given by h) zero forcing 
equalization of the form x e ^ = h(k/T u )~ 1 Xki (or alternatively 
MMSE equalization if the noise variance is known) is per- 
formed where h is the transfer function of the channel (Fourier 
transform of h). 

B. Gabor Multicarrier signaling 

The standard OFDM setup presented in the previous section 
can be embedded into a generalization of MC signaling as 
already proposed by several authors. Because our analysis is 
based on this more abstract formulation, we will give a brief 
introduction. Furthermore we include some remarks on facts 
from Gabor theory that are important for our investigations. 
Hence we focus on a multicarrier system where the transmitter 
modulates data symbols x mn e C on transmitter waveforms 
-f mn with 

lmn(t) = ( S nT ,mF 7 )(*) = l{t - flT)e^ mFt (2) 

being time- and frequency-shifted versions of one transmit 
prototype pulse 7 G £2^). Those sets of functions, denoted 
here as Gabor sets, are generated by unitary representations 
on £ 2 (M) of the so called Weyl-Heisenberg group [11], 
namely the time-frequency-shift operators. Their definition is 
obviously not unique, but (2) is a valid choice. Most of the 
calculations later on can be done by using 

—i2ixab c 

°a,b — e >->-a,-b 

Q Q — i2-nad q 

'~>a,b'Jc,d — C "a+c,Hii 

S a , b S c4 = e- a < ad - bc ^S c , d S a , b 

The rules can be easily verified and are essentially the Weyl- 
Heisenberg group operation. The bandwidth efficiency e (in 
symbols) of the signaling given in (2) is e = (TF)~ l . The 
synthesis of the baseband transmit signal corresponding to the 
transmit symbol sequence x = (. . . , x mn , . . . ) T is performed 
via the Gabor synthesis operator T related to the pulse 7. 
This operator is defined as 

^X ^ ^ ^mn7mn ( ^ ^ X mn S ' nT,mF)^f 
(mn) (mn) 



We will call a point (mn) € X from now on TF-slot. It 
represents in analogy to OFDM the mth subcarrier of the 
nth multicarrier symbol. The index set X <Z 1? itself again 
formally refers to the subset of the rectangular lattice FZ x TZ 
on the time-frequency plane used for transmission. Without 
loss of generality we can embed each x e C 1 into C z by 
setting x mn = for (mn) e 1?\X. Hence the transmit signal 
is given as 

s(t) = (Tx)(t) (3) 

which essentially represents a transmitter side filterbank op- 
eration. The signal at the receiver, after passing through the 
channel 7i and the time-variant distortion S, is 

r(t) - (SnTx)(t) + (Sn)(t) 

A linear MC receiver projects the received signal onto the 
Gabor set {g m n}(mn)ei t0 gi ye trie sequence x. This can be 
formally written by employing the Gabor analysis operator 
G* that corresponds to the pulse g, i.e. 

x = G*r =(..., {g mn ,r), ...) T 

where (x, y) = x*y is here the standard inner product on 
£2(B0 and the operation •* means conjugate transpose. It is 
easy to see that G* is the adjoint of G (same for T* and T). 
Thus, the operator G* implements the receiver side filterbank. 
Now the overall transmission chain is given as follows 

x = G*[SHTx + Sn] = G^SUTx + G^Sn 

H n 

Some important properties can directly be stated in terms of 
synthesis and analysis operators: (1) perfect reconstruction of 
the data symbols is expressed as G*T = I (biorthogonality), 
(2) orthogonal transmitter waveform design means rT = I. 
The matrix T*T is called the Gram-matrix of the transmitter 
pulses. Several properties can also be related to the operator 
IT*. For example in an orthogonal transmitter waveform 
design IT* is an orthogonal projector on the transmitters 
signal space. If the operator norm of IT* respectively TT 

b,= sup ||rvill/||/||^ 

fec 2 (m 

is finite the set {j m n}( m n)ez, 2 is a Bessel sequence and B 1 is 
called the corresponding Bessel bound. If in particular there 
holds < A y — inf /G £ 2 ( R ) ||r*/|j|/|j/||2 the set establishes 
a frame [12] (in our definition it is a frame for £ 2 (R)) and 
IT* is called the frame operator. Furthermore it is a Gabor- or 
Weyl-Heisenberg frame due to the underlying group structure 
(for an introduction see for example [13]). An important case 
arise if A 7 = _B 7 , i.e. IT* = _B 7 L Then this establishes a 
tight frame, which can be understood as a generalization of 
orthogonal bases to overcomplete expansions. 

Finally we adopt the following normalization of the pulses. 
The normalization of g will have no effect on the later used 
system performance measures. The normalization of 7 is 
determined by the transmit power constraint. In our system 
model this can be absorbed into noise scale, thus we assume 
g and 7 to be normalized to one. Furthermore we assume 
a noise power of a 1 per component for the projected noise 
vector h and ~E{xx*} = I. 



III. Interference Analysis 

For this section we provide a rather generic approach for the 
evaluation of the desired performance measures. The reason is, 
that some drawbacks of pure textbook OFDM can be overcome 
with optimized MC transmission schemes. A lot of OFDM 
evolutions have been proposed, where namely pulse shaping 
and different time-frequency densities and constellations are 
considered. Motivated by this observation we start with a 
generic theorem, followed by a corollary related to our channel 
model. Both handle a large class of linear distortions in 
multicarrier transmission. Only in the last step we restrict 
ourselves to OFDM. 

General interference analysis: Writing the received complex 
symbol x k i in the absence of AWGN yields 



xm —HkiXki + (HkiM — Hki)xki + 

V Hu x ( 5 ) 

/ J 11 kl : mn^mn 



H/(M) 



ICI 

def 

where we defined the following expectation value Hki = 
^s{Hki,kl}- Further we define the second moment with 
respect to the statistics of S as Pki = ~E>s{\Hki.ki\ 2 }- Thus 
the transmitted symbol Xki will be multiplied by a constant 
and disturbed by two zero mean random variables (RV). The 
first RV Aki represents a distortion which comes from the 
randomness of S. This part can be understood as a noise 
contribution if the receiver does not know H k ^ki- But with 
proper tracking of Hki.ki tne receiver can "move" the power 
of A kl , given as D 



ki = Pki - \H k i\ 2 , to the desired signal 
contribution yielding an improved performance. Examples are 
the correction of the common phase errors and Wiener phase 
noise tracking. For many applications tracking is mandatory 
to ensure an overall system performance, where interference 
cancellation has less priority due to the steeply increasing 
complexity. Therefore the RV ICI (interference from other TF- 
slots, thus intercarrier and intersymbol interference) remains 
and gives a noisy contribution of power I kl . 

Using this argumentation we establish two performance 
measures important for uncoded communication. With the 
separation in (5) we define the signal-to-interference-and- 
noise-ratio in the TF-slot (kl) for the case where Hki,ki is 
exactly known at the receiver - namely SINR^. This measure 
represents the signal quality if the receiver performs ideal 
tracking of the single TF-slots. If performing no tracking we 
will represent this with a remaining sinr^ < SINR^. These 
definitions and the corresponding bounds are summarized in 
the following theorem. 

Theorem 1 If the realizations HkiM are known to the re- 
ceiver, the SINRfcz of the TF-slot (kl) is lower bounded by 



SINR fe; 



def P k l 



> 



Pkl 



o 2 + hi ~ ° 2 + B~0u - P, 



kl 



(6) 



If it is only possible to track the mean Hki, the remaining 
sinr^ < SINRfc; is lower bounded by 



def 

sinr fe( = ■ 



\H 



kl 



D, 



> 



kl 



a 2 + B 7 p H - \H 



(7) 



kl 



where (3 k i = Es{\\7~t* S* gkiW^} denotes here the channel 
bound and B 7 is the Bessel-bound of{-f mn }. Equality is given 
if the set {jmn} is a tight frame. 

Essentially we used the principle of energy (power) conser- 
vation here to upper bound the received power in each TF-slot. 

Proof: Let 1 C Z 2 . First observe the following upper 
bound 



E 5 { 



E 

(mn)eZ 2 



Hu,mn\ 2 } = E s {(g ki , SH.rr*H* S* gki)} 



< B, ( E s {\\S*H*g k i\\ 2 2 } = B^u 



(8) 



where £? 7 is the Bessel bound of the Gabor set {7 mn }(mri)ez 2 ^ 
i.e. IT* < Bjl Equality in (8) is achieved if J = Z 2 and 
the Gabor set establishes a tight frame. Inequality occurs first 
due to the possible incompleteness of the set {7/d}(fc()ez 2 - 
Then the adjoint channel could map g k i into the orthogonal 
complement of span{7 mn }(„ m ) eZ 2. Moreover (8) can be seen 
as "uniformity" property (related to the condition number) of 
the mapping IT* . Furthermore there could be a sampling loss 
ifJ^Z 2 . 

We can now write the total amount of received signal power 
in the TF-slot (kl) as 



E x , n ,s{\xki\ 2 } 



E 



^s{\Hkl,mn\ 2 } + & 



(8) 

< B 1 (3ki 



(9) 



\R k i(SH)\ 2 + a 2 < B^i + a 2 

The second inequality is caused by the limitation to X C Z 2 . 
The associated contribution, given by 



\Rkl(SH.)\ 2 = E 5 { ^2 \Hkl,mn\ 2 } 1 = Z (JO) 

is the interference from non-existing TF-slots that were 
counted in the infinite sum in (8). Its computation would 
improve (9) in particular at the boundaries of 1 but in the 
following we will neglect this term. Finally, based only on 
the "gain" P k i and the channel bound 0ki, we arrive at the 
following bound to the ICI-power I k i 



hi < B 7 f3 k i - P] 



ki 



(11) 



which then straightforward leads to the bounds on SINRfei and 
sinr fei . ■ 
The latter theorem is a selection of the worst case scenario 
where all symbol energy sent by the transmitter is uniformly 
collected at the receiver. The calculation of the lower bounds 
is in most cases simpler than a direct study of the interferer. 
The importance of Theorem 1 relies in the fact that it is of 
very general type in the sense that concepts like (bi-) orthog- 
onality and completeness of neither the transmit sequences, 
the receiver sequences nor jointly are required. So it is well 
suited for studying distortions that can not be formulated 



within orthogonality of the subcarriers. Moreover, it provides 
a tool for performance evaluations for general non-orthogonal 
multicarrier schemes. 

Using these bounds requires the computation of the Bessel 
bound B Jf which is independent of Ti and S and only 
related to the fixed transmitter setup. For example if {^ m n} 
is an orthogonal basis for its span it follows that IT* is 
the orthogonal projector on span{7„ m }( mn ) eZ 2, i.e. B 1 = 1 
is the minimal achievable Bessel bound for {"f mn } being 
all normalized. For overcomplete sets the minimal Bessel 
bound, achieved by tight frames, is given by the redundancy 
introduced by the normalized {^ mn }. For Gabor sets we have 



B- t 



max{l,— } 



(12) 



Given _B 7 , it remains to compute Hki (or Pm) and 0m for 
the Ti,S. By observing that B 1 is related to the transmitter 
only, we have the desired separation between the system setup 
and distortion. Note that distortions like frequency offset and 
phase noise are obviously not completely known neither at the 
receiver nor the transmitter, so in practice we have to study 
sinrfc; given by (7), hence assuming that at least Hm is known. 
This implies that the receiver corrects the phase so that the 
"full power" \Hki\ 2 can be used for signal reception. When 
we perform an average over the channel we indirectly use 
(6) because we have to assume the channel must be ideally 
known for equalization. Moreover the asymptotic performance 
for ideal tracking of S can be obtained from (6), as shown later 
on for Wiener phase noise. An explicit use of (6) is given in 
[14]. Finally, if the distortion S is not random, it is SINR fc; = 
sinr fe( . 

Incorporating the time-invariant channel: In the aim of 

using (7) of Theorem 1 for the case where Ti is time- 
invariant, known to the receiver and S represents the time- 
variant distortion, we have to specify f3 k i and Hki- Moreover 
it is straightforward to extent Theorem 1 and carry out the 
average over the channel. Hence, let Ti given as 



Ti= [ h{ T )S T ,odT. 
Jo 



(13) 



where h is a realization of the (causal) channel impulse re- 
sponse h with finite maximum delay spread t^. A common sta- 
tistical model for a time-invariant channel is E,{h(ti)h(t 2 )} = 
Ph(ti)S(ti—t2) where p/j is the power delay profile and \\ph\\i 
is the overall channel power (path loss). For this scenario an 
evaluation of (3m (Appendix B) yields 



/3 H ^ f E s {||^ < S* 5H |l2}<l|/i||^E s {|| < S^||^} 
<r d \\h\\ 2 2 Es{\\S*g kl \\ 2 2 }=T d \\h\\l 

and averaging over the channel 

V-h{M = v-H,s{\\n*s*g kl \\ 2 2 } 

= n Ph (T)Vs{\\S* Tfi S*g kl \\ 2 }dT 
Jo 

= ||Pfc||iE s {||5* flfc ,||l} = ||pfc||i 



(14) 



which is independent of (kl). Here we have assumed that 
E 5 {||5/||2} = || /|| 1 which covers the frequency and timing 
offset as well as phase noise. The effective channel matrix for 
a fixed realization h, given as 



(13) 



h(T)(g kl ,SS T , j n 



dr 



(15) 



could still be a RV. To separate the channel from the distortion 
in the evaluation of H k i as much as possible let us define 



Es{(5fcz,SS T ,o7fci}} 



-i2ix k FT 



-i2TTkFT 



Es{(ff> Sit^fSSit^fStaI)} 
(g, ^s{S^ T , kF SSiT, k FS Tfi }'y) 



(16) 



= e- i2 ^ Skl (r) = (S* 0[ <kF s kl )(r) 



which essentially contains the distortion of the rth path 
contribution in terms of the pulses conjugated by Sir, k F, i-e- 
"shifted" to (IT, kF) in the time-frequency plane. The mean 
diagonal (with respect to S) is then given as 



H 



kl 



h(T)F,s{(gki,SS Tfi 'y k i)}dT 



(So.kFh, Ski) 



(17) 



For the channel average we need the second moment of (17) 
with respect to Ti as already intended in Theorem 1 because 
the channel is known, thus 

V-h{\Hm\ 2 } - En{\E s {(g k i,STi lk i)}\ 2 } 

^(T)|E 5 {(.9 fe; ,55 r ,o7fc()}| 2 ^ 



def 



Ph{T)\s k i[r)\ 2 dr = (p h , \s k i\ 2 } 



The disturbances are bounded now as 

hi + Dm < B 7 (3 k i - \ (So, k Fh, s k i}\ 2 
En{hi + D M } < B 7 ||p h ||i - (p h , \s kl \ 2 ) 



(18) 



Let us summarize the refinement of Theorem 1 in the follow- 
ing corollary 

Corollary 1 If the channel realizations of Ti are given as a 
convolution with the impulse response h known to the receiver 
and the distortion S is only known in the mean, the sim k i is 
lower bounded as 



H 



kl 



Smikl <j 2 + Im > <? 2 + B 1 j3 kl - \ (S ^ kF h,s k i}\ 2 ' 



(So, k Fh, Ski) 



(19) 



Furthermore in the average over the channel statistics it 
becomes 



V-h{Hm} 
sinr fc; = — : ^ > 



{Ph, \sm\ 



a 2 + -E H {I kl } ~ v 2 + Bjp h \\ 1 -(p h ,\s kl \ 2 )' 



(20) 



where ph is the power delay profile of the channel. 



So it remains to compute sm(t) (independent of h) and j3 k i 
(depending on h) to get a "worst-case" sinr fer bound. We will 
mainly concentrate on the channel average where it remains to 
compute \sm(t)\ 2 only. The phase of Sm(t) is also important 



to get a view of what the receiver has to correct - separately 
or within the channel equalization. 

cp-OFDM specifics: Before proceeding by applying Corollary 
1 to the problem of time-frequency offsets and phase noise, 
we will introduce a slight modification for cp-OFDM. The 
OFDM transmitter does not exploit an orthogonal set when 
using a cyclic prefix. In the Appendix A it is shown that 
the Gram matrix r*r is block-Toeplitz with the maximal 
eigenvalue given as twice the bandwidth efficiency e, thus 
B 7 = 2e. By using this value Corollary 1 covers an arbitrary 
linear distortion. However it can be shown, that for distortion 
considered in the paper eB g instead of B 7 can be used, where 
B g = 1 is the Bessel bound of the orthogonal receiver set. 
The latter is related to the cyclic structure inserted by the 
special choice of g and 7 in cp-OFDM. This improves the 
prediction in (11), hence we define for our application now 

IV. Timing and Frequency Offsets 

Performance evaluation of communication systems under 
timing and carrier frequency offset is of fundamental impor- 
tance. In particular OFDM systems suffer from a mismatch 
of local oscillator frequency at the receiver with respect to 
the carrier frequency. That means that the decoupling of 
the subcarriers in time-invariant channels achieved with the 
cyclic prefix OFDM is destroyed. First let us start to establish 
formulas for the general offset problem and then refine them 
to OFDM. 

A. The General Offset Problem 

As a simple application of section III we will now study 
synchronization errors that consist of a non-random time- 
frequency shift, i.e. S = Sd,v Because S is non-random and 
unitary it follows that ||<S*gfcz||2 = Hsjwlh = 1. We have to 
evaluate (16), thus 

s kl (r) = e i2 ^ lT - dkF Hg,S T+d ,^) 



(21) 



which is related to the cross ambiguity function A S7 (-, •) of 
the pulses g and 7 



A 97 (d,i/) = (g,S d , v ')) 



(22) 



It will be helpful later on to introduce the normalized offsets 

def ~ def 

v = v/F = evT and d = d/T = edF, such that we have 

Skiir) =e 427r ^-^/ £ A 37 (d + T, v) 

Before proceeding further, we like to state that in the absence 
of a channel we already have 

Hm = Sfcj(O) = e i2 ^ l - ik ^A gi (d,u) 
(3m = 1 (23) 
hi < B 7 - |A 97 (cf»| 2 
for arbitrary d and v. Using this in (7) yields 

|A 37 (d, v)\ 2 



sinr > 



a 2 + B 7 - \A g7 (d, v)\ 2 



(24) 



which is achievable if the receiver ideally corrects the common 
phase error, hence the phase of Hm- Then it turns out that 
SINR fe( = sinr/ci since this distortion is non-random (D k i = 0). 

Let us consider the case where a channel is present. Com- 
puting 



\{s ,kFh,s k i)\ 2 

(Ph, \s k i\ 2 ) 



/ Hr)e 



-i2nTkF 



A gl (d + t, v)dr\ 



Ph{T)\A gl {T + d, v)f 



(25) 

and using Corollary 1 gives first the bound on sinr^ for a fixed 
channel. The second equation gives the bound on sinr in the 
channel average, which turns out to be independent of (kl). 
Both hold for general (Gabor-based) MC schemes given by 
its pulses g and 7. To establish an improved channel bound 
Pki with respect to the pure bound already given in (14) is 
quite difficult for this general constellation. However for the 
frequency offset alone we get (given in Appendix C) 

Pki < \\H* 9kl \\ 2 + 4*^^1(1 + ^1*) (26) 

providing a separation of the frequency offset and the channel. 
Another approximation (not a strict bound) was already pre- 
sented in [15]. In the following we will refine to the problem 
of the frequency offset in OFDM which is of more practical 
importance. 



B. Frequency Offset in OFDM 

Let us consider now a cyclic prefix based OFDM trans- 
mission (instead of B 1 we use now _B fd m = e) distorted 
by constant unknown offsets v and d under ideal channel 
knowledge. The cross ambiguity function for 7 and g as 
introduced in (22) can be compactly written by 



Jcp 



Jcp 



T < 
< T < T, 



cp 



T, 



Cp 



t>T, 



cp 



as 



A gl (T,v) =Ve 



sin7T v{T u 



Jcp I 



(27) 



(28) 



X[-T„,T„]([r]cp)- 



The phase <j)Q = irvT u is related to our choice of time origin 
to = in (1). 

The signal quality in the presence of time- and frequency 
shifts can now be directly obtained from (28). Apart from 
[ ] cp and yj~k (the loss in mean signal amplitude due to the 
cyclic prefix) (28) agrees with the well known auto ambiguity 
function for rectangular pulses g = 7 of width T u . If the 
system exhibits a time offset d only and [d + Td] cp = 0, the 
time dependency in the cross ambiguity function cancels, thus 

Ski (r) = e -^ ik l f - A gl (d, 0) = ^ e <[*>-2*&/e] 

and only phase rotations occur (normally corrected by channel 
estimation and equalization). Contrary to this, time offsets 
with [d + T d ] cp 7^ causes interference. For frequency offsets 
interference occurs immediately as seen from Fig.l. Going 



cross ambiguity function A nv (T,v) tor cp-OFDM 
SNR loss 




: [samfsfes] 



v [samples] 



Fig. 1 . Ambiguity function for cyclic prefix OFDM - The ambiguity function 
A s7 (t, u) describes the behavior of the pulse shaped system with respect to 
single time-frequency shifts, hence natural arise in the offset problematic. It 
illustrates the cp-OFDM fundamentals that the magnitude stays constant for 
time offsets r with [r]cp = 0, where it is rapidly decreasing in v yielding an 
increased interference. 



back to the case [d + Td] cp = 



/i^>o k k/T u ) 



H H = (S , kF h,s kl ) = e «*(«-*0A 

(29) 

holds. Obviously v (d) induces a rotating phase over the time 
slots I (frequency slots k) as one would expect. With (20), 
(26) and (29) we get in the channel average 



> 



sin nv 



/{-Kvf 



o 2 /{A\Ph\\i) + {l-sm^v/{™y) 



where we 

\A gl (T + d,v)f 



again used the 

= |A 97 (0,^)| 2 for [r d + d] c 



fact 



that 
0. This 



result is consistent with a lower bound presented in [3]. 
Restricting to P < i (half the subcarrier spacing) as has 
been done in [3] we immediately obtain their result together 
with an analytical expression for their numerically estimated 
bound on the interference. It can be found by observing that 
INT < (1 - 4f) sin 2 irv = 0, 5947 sin 2 ni), thus 



sinr > (30) > 



1 + 0,5947 sin 2 ttv ■ WphWi/o 2 



(31) 



Note that this bound holds only for v < 0.5 and is less tight 
than (30) (see Fig.2). 

V. Phase Noise 

The continuous time system model we consider in this 
section is 

r(t) = (SHs)(t) + (Sn){t) = d(t)(Hs)(t) + 6{t)n(t). 

Obviously the multiplication by the phase noise process 9{t) = 
e i4(t) fulfills E s { US* /|| 2} = || /|| 2 and does not change the 
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Fig. 2. SINR due to Frequency Offset in OFDM systems - The lower 
bound on the SINR = sinr is shown over the normalized frequency offset 
v = u/T u and compared to the lower bound obtained in [3]. The latter is 
a numerically fit which is only correct for u < 0.5. The bound presented in 
this paper posses the correct behavior SINR — > as u — > 1. 



noise statistics. To derive sinrj.; we start with (17), respectively 
(16), which gives 

su{t) = (g, Si TkF eSiT,kFS Tfi j) 

= (9, S* T0 6SiT,oS T fij) (32) 

= {9, S-lTfiOSlTfiSTfij) 

where the first moment is defined as 9{t) = {#(£)}. The 
second step follows because 9 is a pointwise multiplication, 
so that the frequency shifts will cancel out. 



(30) A- Gaussian Phase Noise 



A typical model which occurs in phase synchronization 

loops is 4>(t) - 7V(0, S0) with E{>(t)0(i + r)} = C^(r) = 

/ S^(f)e l27T ^ T df. It was already observed in [2] that sinr/^ is 

independent of the phase noise spectrum S<f,(f), where the 

authors only considered an classical OFDM system based 

on rectangular pulses (without cyclic prefix and additional 

channel). For the general bounds presented in this work this 

is a direct consequence from the fact that the bounds depend 

— _ 

only on the first moments. Thus, with the mean 9{t) = e 2 
for the Gaussian case we get 

Skl(r) = e~~ (g,S T fij) = e~~A ff7 (r,0) 

and we continue (except of the constant factor) as in Sec. IV. 
That is 



Ha =e 



h(r)e 



-i2irkFT 



A 97 (T,0)(ir 



If no channel is present we obtain for general (Gabor-based) 
MC schemes with A 97 (0, 0) = (3,7) 

. > I (9, 7) I 2 

Smr " eM^ + ^-I^T)! 2 
and for OFDM (using £? fdm = e) 

1 

sinr > 



e^(o- 2 /e+ 1) - 1 



For the channel average we get similar to (25) 

{Ph,\ski\ 2 ) =e~ s * f d Ph (T)\A g7 (T,0)\ 2 dT 
Jo 

determining sinr (see Corollary 1) for general (Gabor-based) 
MC schemes. For OFDM transmission (using B f dm = e and 
(28)) this gives 



\h(k/T u ) 



e s*(aye + p kl )-\h(k/T u )\2 



if the channel delay spread does not exceeds the cyclic prefix. 
While in the latter the calculation of flki is still left open 
for practical applications, we obtain in the channel average 
directly 

1 

Smr " e s ^ 2 /MMi) + l)-l 



B. Wiener Phase Noise 

A widely used model in frequency synchronization is 



4>(t)= j $ 



tfidr 



with the instantaneous frequency j^4>(t). The power density 
spectrum (PDS) of the signal r(t) corrupted by phase noise is 
given by 

S r (u>) = (So * S-hs){uj) 

with S e {u) = J^e-^^Ceir) and C e being the auto- 
correlation of the phase noise process. With the application 
of the Wiener-Khintchine-theorem the autocorrelation can be 
expressed using S^ (the PDS of the instantaneous frequencies) 

C (t) = e~£ ^ s «(") sinc2 (tt) 

For S^(w) = = const the autocorrelation of the process is 
given by 

<7,(r)=e-^M 
with the typical Lorentzian PDS 

Se(u 



Si + 4(2trj) 2 



In the presence of Wiener phase noise communication via 
coherent detection is not possible due to the infinite distortions 
of the phase. A common approach is to correct the phase from 
period to period, therefore we use 

9{t) = e^*-**-) 

with </>(•) being a realization of the Wiener process <p(t) = 
j Q (j>dr on [0, oo) and t sync denotes the time of the last phase 
synchronization. For simplicity let us assume that i sync = 
kyncT, i.e. a multiple of the symbol time shift. The mean of 
9(t) (defined on [i sync ,oo)) is 

9(t) = e-^(*-**-) = e ^H(i) 
where 0\(t) is a phase noise process defined on [0, oo). 



We will need the time-frequency shifted version 
of 6\ in (32). Thus with the commutation relation 

S-iTfi~6i = e^ lT 9i S-it,o this gives 

sm{t) = (g, S-iT,o6SiT,oS Tt oj) 



e" ts " c (s, S^TfiOiSiTflSr.a-f) 
e¥(^-')( 9 ,5 1 S Ti07 > 



(33) 



which depends obviously on the pulse shapes g and 7. The 
last step is correct in a rough sense only. The reason is that 
the pulse shapes could be much longer than T. So one has 
to assure that the domain of 9\ is not/or marginally violated 
(this is obviously not relevant for cyclic prefix OFDM). We fix 
the time origin now so that Z sync = and normalize the phase 

def 

noise power with respect to the subcarrier spacing p = S^/F, 
i.e. 



— l {g,0 1 S Tfil ) = e-^{g 1 9 l S Tfil ) 



Now, using (33) we can directly establish the following: if no 
channel is present a bound on sinrfe; is 



sinr fei 



> 



I {9, on) I 2 



ePi/^+B^-lig^nW 
and in the channel average 

Io d Ph(r)\ (g, OiSrfil) \ 2 dr 



(34) 



sinr > 



e^A( a 2 + B 7 || Pfc || 1 ) - ^ p^g^Sr^dr 

(35) 

Both formulas hold again for general (Gabor-based) MC 
schemes. And - from (34) and (35) one can see that there 
is an inherent exponential degradation of the signal quality. 
To obtain closed formulas it is left to calculate (g,9iS T .oj)- 
For OFDM we get 



Ma-e-^( r --l^l)) X[ _ T ., r j([r],) 
With the normalized phase noise power p = SiT u 



(5,0iS T ,o7> = i^O 

>5 1 1 u 



<<7ASt,o7> = — 
P 



e"«) 



follows as long as the channel delay spread does not exceed 
the cyclic prefix. If no channel is present (34) yields 

1 

4( l- e ->/»)» (' 7 / t + 1 ) _ 1 

The bound for the channel average given by (35) reads now 

1 

sim ~ P 2 e»'/' — / 21/ 11 11 \ -< \ — 7 (37) 

4(i-^/')» C* 7 /( £ lblli) + 1 )- 1 

And finally for a particular channel realization h we obtain 

\h(k/T u )\ 2 



sim k i > 



{ £:%> (<r 2 /e + M-\h(k/T u )\< 



(38) 



In the last estimates we used again that i? fd m = e - Note 
the exponential decay in (38), thus tracking for Wiener phase 
noise is crucial. One can directly obtain the tracking gain, 



hence I = (I = 1) means that phase synchronization 
was performed in the current (previous) OFDM symbol. It 
is interesting to see which asymptotic performance results if 
the receiver ideally removes the phase noise on each subcarrier 
but does no further interference cancellation. We can give an 
answer to this question by using formula (6) from Theorem 1, 
namely the bound on SINR^;, which is shown in the Appendix 
D to be 

1 

SINRw > - 2 — ; > sinr fe; 



r (a 2 /e+l)-l 



The graphical summary 
(36), is shown in Fig. 3. 



4(p-2-2e-P/ 2 ) 2 

i.e. the comparison of the latter to 



SINR 




Fig. 3. Wiener phase noise tracking in OFDM systems - The sinr; for 
I = 0,1,2 is shown over the normalized Wiener phase noise power p 
and compared to the asymptotic performance given by SINR. The latter 
corresponds to a perfect tracking of phase noise in each subcarrier, but without 
a cancellation of the induced interference. 
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Fig. 4. MSE due to Frequency Offset in OFDM systems - the impact of the 
frequency offset and its prediction on the MSE over the normalized offset 
with and without a LTI channel. 



large enough so that the main degradation is i/y with many 
contributions (central limit theorem). 
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VI. Performance Evaluations for OFDM 

A. Frequency Offset 

To demonstrate the prediction of the degradation due to 
a frequency offset we present in Fig. 4 the theoretical and 
simulated symbol estimation error MSE = f \xfl — Xki\ 2 over 
the normalized frequency offset for a fixed sample channel. 
With the approximation in (26) for the channel bound flki the 
theoretical prediction agrees nicely with our simulation. An 
unknown frequency offset has significant impact on the MSE 
especially if a (known and equalized) channel is present. 

B. Phase Noise 

As an example of the evaluation in section V-B we present 
here the symbol error rate (SER) of a cyclic prefix OFDM 
system in the presence of Wiener phase noise. For the SER- 
prediction the interference due to phase noise is assumed to 
be Gaussian (see Fig. 5). Then it can be treated as additional 
noise. We like to state that this is only appropriate for S) 



Fig. 5. BPSK symbol error rate - the simulated impact of Wiener phase 
noise on the BPSK performance and its prediction is shown. The phase is 
synchronized at each OFDM symbol. 

VII. Conclusions 

We derived a framework for the evaluation of bounds 
on uncoded system performance of linear distorted (Gabor- 
based) multicarrier schemes with inclusion of a time-invariant 
channel. We identified the dominating terms determining in- 
terference levels for the case where the receiver has perfect 
knowledge or only mean knowledge on the distortion. Our 
contribution provides analytical insights into the interrelation 
of pulse shapes and time-frequency density of Gabor systems 
to practical problems of imperfect, hence distorted radio fron- 
tends. The bounds apply without requiring (bi-)orthogonality 
of the subcarriers commonly needed for those evaluations. Our 
study was motivated by the impact of time-variant distortions 
on the OFDM performance caused by imperfect receiver 
structures. Therefore we applied the theoretical framework 



on time-frequency offsets and phase noise, both being effects 
limiting the performance of current OFDM implementations. 
Finally we verified our theoretical predictions with computer 
simulations. 



Appendix 
A. Bessel bound for cyclic prefix OFDM 

The Bessel bound is given by the norm of IT* which 
is equal to the largest eigenvalue of the Gram matrix TT. 
Computing this for j(t) = -^= X [- Tcp ,T u ](t) is 



(r*r) fcI , 



Si n e 



sin ^ (to — k) 
f (to - k) 



where e = T U /(T U + T cp ). This is a Toeplitz matrix in the 
frequency slots k and m generated by the symbol 



iir(2u!-±)n 



sin — n 



— 71 



2e ^ cos7r(2w — - )n ■ sin -n 
= 1 + — > — 

n=l 



1 



= 1 + e- e[(l/e-uj) mod 1 + uj] = e(|_- - wj + 1) 

(39) 

where in the last step we restrict ui £ (0, 1). Its known that the 
spectrum of the infinite Toeplitz operator T*T is dense in the 
image of </>. Thus the Bessel bound is given as £? 7 = H^Hoo- 
For cyclic Toeplitz matrices already for finite dimension N the 
kth eigenvalue is given as <p(k/N). As seen from (39) </>(u) 
is a step-like function taking on [0, 1] only the two values 
e( [ij + 1) and e(|_| - lj + 1). A special case is e = 0.5 
where the first value is not anymore in (0, 1). In that case 
and also for e = 1 the spectrum is constant <p(u>) = 1, hence 
the set {"f mn \ forms an orthonormal basis for its span. For 
e > 0.5, relevant for the application in cp-OFDM, <f>{uj) is 
a step-function with the levels e and 2e, thus in general we 
have to use B 1 = 2e. Note here that normally (but not in this 
paper) the transmit pulse is normalized for transmit power 1, 
thus ||7|| 2 = l/e and then follows B 1 = 2. Using B 1 = 2e 
in (18) gives the most general bound on the disturbances. In 
particular this is needed if S represents non-causal operations, 
like S- T . v with t > 0, occurring in time-offset correction. 
For the distortions considered in this paper one can show, 
that one can equivalently use eB g where B g is the Bessel 
bound of the receiver set {g m n}(mn)ei, 2 - Its easy to verify that 
(mn)£i 2 is an orthonormal set, i.e. GG* an orthogonal 
projector onto its span, hence B g = 1. 



B. General Channel bounds 

n* is a convolution with h R {t) =h<~^). With (h R y = h 
this gives for the channel bound /3ki 

Pti - E s {\\(S*g kl ) * h R \\ 2 2 } = E s {\\h- {S*g kl y\\l} 

<m/ i | 2 ii P .E 5 {|ii(5*. 9fe; ri 2 iij 



The last step is for 1 = l/p + 1/q (Holder's inequality). From 
practical point of view the essential support of h is often less 
then g. Therefore using p = 00 and q = 1 



^<||/ i || 2 00 -E s {|| < S^|| 2 } = ||^|| 2 <) 
with E 5 {||5*/||} = 11/11 and \\g kl \\ 2 = 1. 

C. Channel bound for the frequency offset 
Let S = So,„ and [i = h R , then 

V s {\\n*S*g k i\\ 2 2 } = J \»(u)\ 2 \g kl (u;-v)\ 2 du; 

= J + v)\ 2 \g k i{u)\ 2 <kj 
We can express the first part of the integrand as 

|A(^)| 2 = |AMI 2 + | r + %'(/)d/l 2 



(40) 



(a) 



+ 25R{A(c) f U+V p!{f)df} 



(b) 



By using 



\H'(f)\ 2 = \§- f £ d dMr)e i2 ^\ 

= \2tt f * dTTn(T)e l27TfT 
Jo 



< {2*) 2 T d 2 \Kf)\ 2 

we upper bound term (a) as (using Jensen's integral inequality 
and the measure df/\v\) 



(a) 



=M a ir + "?A'(/)i 



- H U R IA ' (/)|2 

Further for the term (b) follows 

(b)<2|£(w) / + dffi'(f)\ 



<2\fi(u)\ f + ' 'rf/|A'(/)| <47r||A|UIAM|r d |H 
< 4 7 r||A|| 2 .r d |H 

Putting both together and using that HAUL < T d||Mll2 gi ves 

|A(w + v)\ 2 < IAMI 2 + 47rr d |^|||A|| 2 (l + T d \v\*) 
< |AM| 2 +47rT d 2 |^|(l + r d |^K)|| M ||| 



D. Wiener Phase noise - tracking asymptotic 

Let S = 6(t) be the pointwise multiplication with the 
phase noise process and Co (r) the phase noise autocorrelation. 
Moreover, for simplification no additional channel should be 
present, i.e. H, = I. Let assume further that the receiver can 
ideally track and correct the phase noise on each subcarrier, 
but does no interference cancellation. The performance of this 
asymptotic situation can be obtained from Theorem 1. The 
SINRfc; in formula (6) is determined only by P^. The channel 
bound is /3ki = 1 (because H. = I) and £> fd m = e (Appendix 
A). Then 

Ph = ~&s{\Hki,ki\ 2 } 

= J J f(ti)C e (ti - t 2 )7(t2)d 2 t 
= J \f(u)\ 2 Se(Lj)du 

Hpf 

with f(t) = g(t)-y(t). For OFDM follows f(t) = 
^X[o,T u ](t), hence |/>)| 2 - e s \f™J? . With the Lorenzian 

PDS Sg(uj) = sT+4(2nuPp ^ or Wiener phase noise follows 



kl 



sin' 



itujT u 4S 1 ; 



(ttu;T u ) 2 S?+4(27ro;) 2 



duo 



4e 



r(5;T„-2 + 2e- — ) 



= ^(p-2 + 2e-S) 

where p — S^T U . Due to "permanent" ideal tracking the (kl) 
dependence disappeared. The SINR bound is now 



SINR > 



1 



4( P -2-2e-P/2)2 ( g2 A+l) - 1 
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